Brane Worlds and the Cosmic Coincidence Problem 
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Brane world models with 'large' extra dimensions with radii in the n ~ 10 — 100 /im range and 
smaller ones at r s < (ITeV) -1 have the potential to solve the cosmic coincidence problem, i.e. the 
apparently fortuitous equality between dark matter and dark energy components today. The main 
ingredient is the assumption of a stabilization mechanism fixing the total volume of the compact 
submanifold, but allowing for shape deformations. The latter are associated with phenomenologically 
safe ultra-light scalar fields. Bulk fields Casimir energy naturally plays the role of dark energy, which 
decreases in time because of expanding r;. Stable Kaluza Klein states may play the role of dark 
matter with increasing, 0(l/r s ), mass. The cosmological equations exhibit attractor solutions in 
which the global equation of state is negative, the ratio between dark energy and dark matter is 
constant and the observed value of the ratio is obtained for two large extra dimensions. 
Experimental searches of large extra dimensions should take into account that, due to the strong 
coupling between dark matter and radii dynamics, the size of the large extra dimensions inside the 
galactic halo may be smaller than the average value. 
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I. INTRODUCTION 

The energy content of the Universe is nowadays, more 
than ever, the subject of intense research. The picture 
emerging is the following [1]; the total energy density of 
the Universe equals the critical value p c ~ (2-10 -3 eV) 4 ~ 
(100 jLtm) -4 (i.e. ft ~ 1), non-relativistic dark mat- 
ter (DM) accounts for at most forty percent of it ( 
Q.DM — 0.2 — 0.4), and the rest is in the form of a smooth 
- i.e. non-clusterized - component, named Dark En- 
ergy (DE) (plus a small amount of baryons, S1& ~ 0.05). 
Luminosity-redshift measurements of supernovae la [2], 
as well as the resolution of the so called age problem [3] , 
agree with this picture and point towards a DE with neg- 
ative equation of state W — p/ p, where p and p are pres- 
sure and energy density, respectively. A viable DE candi- 
date is then Einstein's cosmological constant (W = — 1), 
but more general fluids with W < have been consid- 
ered, emerging, for instance, from the dynamics of scalar 
fields with appropriate potentials. The old cosmologi- 
cal constant problem [4] now becomes double-faced. On 
one hand we still have the old puzzle of the cosmologi- 
cal constant - or the scalar field potential- being so tiny 
compared to any common particle physics scale. On the 
other hand - and this is the new face of the problem- 
the presence of DE today in roughly the same amount 
as DM energy density is embarassing. Indeed, the ratio 
between the two energy densities scales as 



Pdm 
Pde 



,3W 



where a is the scale factor of the Universe. So, the ap- 
proximate equality between the two components today 
looks quite an amazing coincidence in the cosmic history 
if DE has If ^ as required by observations. 

In principle, brane world scenarios offer a suggestive 



contact with the DE problem. Indeed, extra spatial di- 
mensions compactified to a size as large as 100/zm have 
been shown to be a viable possibility, provided that no 
Standard Model (SM) field propagates through them [5]. 
It might then be tempting to attribute the observed value 
of the DE energy density to the Casimir energy associated 
to some field propagating in these large extra dimensions 
of size r, i.e. pde — B/r 4 , with B a constant typically 
O(10" 5 - 10- 4 ) (see sect. IV). 

If the radius r is stabilized, then the energy density 
behaves exactly as a cosmological constant, and we have 
made no substantial progress with respect to any of the 
two problems mentioned above. Indeed, we have no clue 
about why the Casimir energy would be the only non- 
vanishing contribution to the cosmological constant and, 
moreover, the equality between the latter and the matter 
energy density would appear as fortuitous as before. 

In order to gain some insight at least on the cosmic 
coincidence problem it seems then necessary to make the 
radius r dynamical, so that it may evolve on a cosmolog- 
ical time-scale. But this turns out to be quite dangerous. 
The 'radion' field, whose expectation value fixes r, cou- 
ples to the trace of the SM energy momentum tensor with 
gravitational strength. Moreover, in order to evolve with 
r/r ~ H (H being the Hubble parameter), it must be 
extremely light, m ~ H ~ 10~ 33 eV. It then behaves 
as a massless Jordan-Brans-Dicke scalar with O(l) cou- 
plings to matter, whereas present bounds are O(10 -3 ) 
[6]. Thus, the geometric explanation of a dynamical DE 
looks quite unlikely (see however [7]). 

The purpose of this paper is to show that the above 
conclusion is indeed premature, and that standard brane 
world scenarios with large extra dimensions may quite 
naturally solve the cosmic coincidence problem. Unfor- 
tunately, we will add nothing new with respect to the first 
problem, that of the vanishing of all present contributions 



1 



to the cosmological constant larger than 0(100 /im)~ 4 . 

The key point in our discussion is the realization that 
besides large extra-dimensions in the 10 — 100 ^.m range 
there could also be smaller ones. Considering for sim- 
plicity only two compact subspaces, each caracterized 
by a single radius, we note that the trace of the four- 
dimensional energy-momentum tensor couples with the 
total volume of the compact space, that is with the com- 
bination r" s r" ! , where r Sy i and n Sj ; are the radii and di- 
mensionalities of the two subvolumes, respectively. If 
we assume a stabilization mechanism that fixes the to- 
tal volume V of the compactificd manifold but not its 
shape, the potentially dangerous combination of radion 
fields associated with volume fluctuations is made harm- 
less, whereas the orthogonal one, associated with shape 
deformation, is decoupled from normal matter and may 
then be ultra- light. As a consequence, n can grow on a 
cosmological time-scale (such that the associated Casimir 
energy B/rf decreases) and at the same time r s shrinks 
so that V keeps a fixed value. 

The coincidence problem is solved if we associate the 
dark matter with some stable state living in the extra 
dimensions, because the resulting cosmological equations 
exhibit an attractor solution in which the Casimir energy 
and pdm are redshifted at the same rate with their ratios 
fixed at a O(l) value. 

If the DM candidate has a mass independent on r s then 
the attractor corresponds to the equation of state of non- 
relativistic matter, i.e. W — 0, which is disfavoured by 
cosmological observations. 

On the other hand, if the mass of the DM particle 
is inversely proportional to some power of r s , as for a 
Kaluza-Klcin (KK) state (todm = mKK ~ then 
a more interesting behavior is obtained. Indeed, on the 
attractor now we have 

Pdm ~ l/^a 3 ) ~ Pde] , 

which, due to the shrinking of r s , decreases slower than 
a~ 3 during the cosmic expansion. Remembering that a 
fluid with equation of state W scales as p <~ a~ 3 ( w+i ) we 
see that the energy densities scale now with a negative 
effective equation of state. Moreover, the values of W, 
ft dm, and OflE depend only on the two parameters, n s 
and 71/ . 

The above results rely on naive dimensional reduction, 
without assuming drastic modifications to the radion ki- 
netic terms as was done for instance in [7]. 

For sake of clarity, we summarize our assumptions 
here. They are; 

1) the vanishing of any contributions to the effec- 
tive 4-dimensional cosmological constant larger than 
O(100 pm)- 4 today; 

2) a stabilization mechanism for the total volume 
(rather than for each single radius separately) of the com- 
pactified space effective at low energies. 

The paper is organized as follows. In sect. II we in- 
troduce our setting and derive the effective action after 



dimensional reduction. In sect. Ill we discuss a mech- 
anism to stabilize the volume of the compact manifold, 
allowing at the same time its shape to vary. In sect. IV 
we will study the dynamics of the modulus associated 
to shape deformations of the compact manifold, showing 
that an (attractor) solution exists on which dark energy 
and dark matter scale at the same ratio and the uni- 
verse has a negative effective equation of state. Finally, 
in sect. V we discuss our results. 



II. THE SETTING 

The starting point is the action 

S = Sbulk + S4+n 3 + S4 ■ 

The bulk action is given by 

R(Gab 



+ A B + . . . 



(1) 



(2) 



where Af = n s + ni and the 4 + Af- dimensional metric is 
given by 

ds 2 = g AB dx A dx B 

= g^dx'dx" + J2 biixf^dy^dy^ , (3) 

where p, v = 0, ...,3, a,,/3j = l,...,m, and bi(x) = 
ri(x)/r®, with r°'s the average values of the radii today. 
As is the bulk cosmological constant, whereas the dots in 
cq. (2) represent extra fields living in the bulk which will 
contribute to the stabilization of the compact volume, as 
we will discuss in the next section. 

Apart from the overall b^s factors, the metrics 
are assumed to be non-dynamical. 

S , 4+ ris is the action for the fields living in all the 4 + 
n s dimensions. In the following, we will only assume 
that there is at least one stable KK state, which we will 
associate to DM. 

The SM fields are confined to our four-dimensional 
brane, with action given by 



S 4 



'sm^^VO + A 4 ] 



(4) 



where V represent SM fields and we have singled out 
the contributions to the cosmological constant from the 
brane, A4. The dimensional reduction of Sbuik is a 
straightforward procedure (sec for instance [8]). We de- 
fine new 'radion' variables, <f> a (x), as ln&j = A ia <p a with 
i = s, I, a = 1, 2, and 



A = 



1 



2 

2+AT 
2 

2+AT 



(5) 



By also rescaling the 4-dimensional metric as 



2 



9fiv — e ' N '^ >1 fjfii' ■ 



(6) 



where CV = \/2N/{2 + TV), we obtain the 'Einstein 
frame' action, 



u. k = / 



1 



167tG 4 



M 2 2 



a=l 



(7) 



where we assumed that both iJ^W) — 0, as in toroidal 
compactifications, and V = G^+n/Gh is the volume of 
the compact manifold today. 

In the same frame, £4 takes the form 

S A = J rf 4 ^e- 2C ^* [C S M{e- C ^~g^) 

+A 4 ] , (8) 

from which we see that only fa, which we might call 
'the volumon', couples to matter. Indeed, from eq. (5) 
one realizes that fa is the scalar field controlling volume 
deformations, since 



The orthogonal combination, fa, controls shape defor- 
mations of the compact manifold and does not couple to 
SM fields directly. On the other hand, it couples to fields 
living in the bulk, or on the 4 + n s -dimensional branc. 
In the 4-dimensional language, fa couples - besides the 
graviton- to non-zero KK states, since their mass terms 
are proportional to 1/rj ~ exp(— A a fa) (i — s,l). 



III. VOLUME STABILIZATION 

The field fa couples to the SM ("IS (\ Jordan-Brans- 
Dicke scalar, i.e. in a purely metric way. As is well 
known, this type of couplings enforces the (weak) equiv- 
alence principle, but is constrained from gravity tests in 
the solar system and in the laboratory [6]. A massless 
fa would modify gravity at an unacceptable level un- 
less CV < O(10~ 3 ), which is clearly ruled out in the 
present model, as CV ranges from 1 (n s = ni = 1) to \pl 
(n s + ni — > 00). It is then necessary to give a mass to fa, 
that is, to stabilize the volume of the compact manifold. 
In order to fulfill the experimental bounds, the fa range 
should be smaller than about a hundred microns, that is, 
m r f >1 <; 10~ 3 eV. At the same time, we would like to keep 
fa massless such that it may play the role of dynamical 
dark energy today, that is, m^, 2 ^ Hq ~ 10 _33 eV. Thus, 
we are looking for a stabilization mechanism which fixes 
only the total volume of the compact manifold, allowing 
the different dimensions to vary accordingly. 

A mechanism with this property has been discussed for 
instance in refs. [9], and is based on the introduction of a 



U(l) gauge monopole configuration in the extra dimen- 
sions, besides the bulk cosmological constant A B and the 
'brane tension', A4, that we have already introduced. 

For definitness, we will first consider a two dimensional 
torus of radii bi and b s , and then we will discuss the exten- 
sion of the mechanism to higher dimensional manifolds. 
The key point is that U(l) gauge fields can have non-zero 
magnetic flux through the two-dimensional surface, and 
the flux is quantized and topological invariant, 



/ 



d 2 n e mn F 



2nk 



(10) 



45 



,54 



where to, n = 4, 5, e 
Adding the gauge kinetic term, 



1, and k is an integer. 



/ 



d^xV^G^F AB F AB 



(11) 



(A, B = 0, ...,5), to S'buik and solving the Euler- 
Lagrange equations for Fab, a solution can be found of 
the form 



1 ILV 



0, 



,B 



b,b 



i<J s 



(9) with d m B = 0. M* in (11) is a mass scale typically 



0(G 6 1//4 ). Inserting this solution in eq. (10) we find that 
flux quantization forces the electromagnetic field to be 
inversely proportional to the 2-dimensional volume, 



B 

so that eq. (11) becomes 

k 2 TT 2 M 2 



kit 
Vhbl ' 



2g 2 V 



_g e -3C 2 4>l 



(12) 



Like the terms proportional to in eq. (7), and to A4 
in eq. (8), eq. (12) is only sensitive to the volume, that is 
to fa, and not to the shape modulus, fa. These contri- 
butions give the effective potential for fa in the Einstein 
frame, which is the appropriate quantity to minimize in 
order to study stabilization*, i.e 

V{fa) = De- C ^ + Ae- 30 ^ 1 + A 4 e- 2C2 ^ , (13) 



with 



k 2 n 2 M 2 
2g 2 V 



D = VAj 



*In the second of refs. [9] they minimized the effective poten- 
tial in the frame with metric g M „, where the field fa appears 
also in front of the Ricci scalar, which should be taken in into 
account in the minimization. As a consequence, their results 
on the mass of the scalar field differ significantly from ours. 
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Notice that A > 0, which ensures that the potential is 
bounded from below. Requiring that the contribution of 
the potential to the cosmological constant vanishes in the 
minimum, we get the fine tuning condition 

A 2 = AAD , 

which forces D > 0, and then Ab > 0. The poten tial has 
a minimum only if A4 < 0, that is, for A4 = — 2^/ AD, in 




On this minimum, the scalar field (f>i has a mass 



< = V"(<pf n )/M^ 



93/2^2 

±-^ (r °M,) 2 (r°M^ 



MV Ml 



(14) 
(15) 



The above mechanism can be generalized to compact 
spaces of higher dimensions. In order to stabilize the vol- 
ume of a A/"-dimensional manifold, one needs a M— 1-form 
C/(l) gauge field, with A/"-form field strength F W> . Then, 
the flux (10) generalizes to the topological invariant 



Mif- 



J 



and the kinetic energy is again inversely proportional to 
the compact volume [9] . Considering again toroidal com- 
pactifications, the ./V-dimensional case gives a scalar mass 



93/2^2 

^-^(r°M*) 2 ™*(rfM*) 2n < 

KIT 



M 2 ' 

p 



(16) 



The bulk cosmological constant is bounded from above by 
the requirement that it does not induce a bulk curvature 
radius smaller than [9], which gives an upper bound 

on rrifa , 



93/2^(2 / M \ 4 ~ D /™i M 4 

1 * tt {rsM *> \Mj Ml ' (17) 



4-6 /m 



m, 



where we have taken k = 1. Assuming r s M+ 1 and 
M* <; lOTcV we see that the requirement m^, 1 ^ 10 -3 eV 
is in no conflict with the bound above for any ni>2. 

Since m < p 1 > H (H being the Hubble parameter) for 
temperatures T ^ M+, in the following we will consider 
the field 4>i fixed during the late time evolution of the 
Universe, and will study the dynamics of the light </>2 
field. When (f>i is constant and massive, the two frames 
related by eq. (6) are equivalent. 



IV. VAMPS FROM EXTRA-DIMENSIONS 

Integrating the non-zero KK modes out leaves 
'Casimir' contributions to the free energy proportional 
to 1/r^; . For instance, a scalar field compactified on 
a circle of radius r with periodic (antiperiodic) bound- 
ary conditions gives a contibution —5.06 • 10~ 5 /r 4 (4.74- 
10~ 5 /r 4 ) [11], whereas a fermion gives 2.02 • 10~ 4 /r 4 
(—1.90 • 10~ 4 /r 4 ). For higher dimensional torii, the nu- 
merical coefficients change, but the order of magnitude 
is the same (on a two-torus with r\ = r2 = r, we find 
-4.87 • 10- 5 /r 4 ( 1.95 • 10~ 4 /r 4 ) for a scalar (fermion) 
with periodic boundary conditions. 

Consistently with our assumption 1) in the Introduc- 
tion, we will neglect the 0(rj 4 ) term, since it must be 
cancelled by the (unknown) mechanism solving the cos- 
mological constant problem. A possible mechanism could 
be supersymmetry in the bulk or, even without super- 
symmetry, a particle content in the 4 + n s brane, can- 
celling exactly the 1/r 4 contribution, in the same spirit 
ofrcf. [12]. 

Thus, we are left with the Casimir contribution from 
the 'large' dimensions, given by 



V ^ = ^= (^4 eXP ' 1 



(18) 



where we have fixed 0i = and B is a O(10~ 4 — 10~ 5 ) 
coefficient depending on the particle content and the di- 
mensionality of the 4 + A/"-dimensional bulk, which we 
require to be positive. 

As we have anticipated, the other ingredient of our 
model is dark matter. We associate it to a stable KK 
state, whose mass scales as 0(1 /r s ). It will be an ex- 
ample of varying mass dark matter particles, or Vamps, 
which were discussed in a different context in refs. [10]. 

The important point is that the cosmological abun- 
dance of this non-rclativistic relic will scale in this case 



as 



Pdm ~r s 1 a 



exp 



n s Af <l 



(19) 



Since the runaway potential, eq. (18), pushes the field 
4>2 to —00, the mass of the dark-matter particle increases 
during the cosmological expansion, and its energy density 
rcdshifts less than for common dark matter. 

The fact that (f>2 is the canonically normalized version 
of a radion field is the reason for the exponential depen- 
dences in (18) and (19). This is of crucial importance for 
what follows, because exponentials, once inserted in the 
cosmological equations, allow scaling solutions, in which 
Pde an d pdm rcdshift at the same rate. 

Defining 



= 4 



n. 



X = 



niAf 

the cosmological equations are, 



ni 



n. 



(20) 



4 



H 



■2 + 'iH(t> 2 = -0V + \p D M 

if Ml ~ 



PDM + -^-4*2 + V 



- a =-^(P DM + 2M ^ 2 - 2V ) ■ < 21 > 



They admit a solution of the form 

-3 



A + /3 



(do = 1) which is such that , 
the equation of state 



W = T^ = - ( 1 + 4 - 

A + (3 V n i 



log a , 

~ PDM ~a^ w+1 \ with 
(22) 



and fixed ratio, 



n 



DE- 



Pto = 3 + M£ + A) 
Pdm + P4> 2 ((3 + A) 2 
ni(3n s Af + 4n s + n;) 
(4n s + n z ) 2 ' 



(23) 



independent on the scale factor a. The solution above is 
an attractor over solution space if 



n s > 



ni(3ni - 4) 
16 - 3ni 



(24) 



Differently from the case of a cosmological constant, or 
even from quintessence models with inverse power law 
potentials, once the attractor is reached VLde and Qdm 
become independent on the cosmological era, thus solving 
the cosmic coincidence problem. Moreover, the equation 
of state (22) is negative, as required by observations. 

The values of W and SlcEon the attractor are func- 
tions of n s and rii only. In Tab. 1 we list the possible val- 
ues of fi_o_E, W and of H to, to being the present age of 
the Universe. We limited the dimensionality of the com- 
pact space according to the theoretical prejudice coming 
from string theory, i.e. W < 6. 



n s 


m 


Qde 


W 


Hota 


I - 5 


l 


< 0.44 






1 


2 


0.83 


-1/3 


1 


2 


2 


0.68 


-0.20 


0.83 


3 


2 


0.60 


-0.14 


0.78 


4 


2 


0.56 


-0.11 


0.75 


1-3 


3 


> 0.92 






1-2 


4 


no attractor 






1 


■5 


no attractor 







Tab.l The values of CIde, W, and H to for different values 
of n s and n; such that Af = n s + ni < 6. 



V. DISCUSSION 

The first noticeable fact about Tab. 1 is that the ob- 
served range for the dark energy density, 0.6 ^ SIde ^ 
0.8 uniquely selects the number of 'large' extra dimen- 
sions to be ni = 2, the same value that is required by the 
totally unrelated issue of solving the hierarchy problem 
with 'millimeter' size extra dimensions [5]. 

Indeed, at the level of Tab. 1, we have not yet inserted 
any information about absolute scales, such as 100 /zm, 
TeV, Ho and so on. We find it a remarkable and inspiring 
'coincidence' that the observed balance between £Idm 
and OflE points to the same value, n; = 2, obtained from 
scale dependent considerations, such as reproducing H 
or solving the hierarchy problem. For n; 2 we either 
find that the attractor corresponds to energy densities 
outside the observed range or that the couplings A and 
(3 lie outside the limit of eq. (24) and correspond to a 
different -unphysical - attractor. 

The values for H to listed in the table are obtained 
from 



Ho(to — tatt) 



3(W- 



where t at t and a a tt are the time and scale factor when 
the attractor is reached, assuming a a tt < 10~ 2 . Taking 
the 95% confidence level lower limit on t from globular 
cluster age estimates, to > HGyr [3] we obtain the lower 
bound 

HoU) > (0.71 - 0.79) 

where we have varied H in the range (63 — 
70) kms _1 Mpc _1 . The above bound is inconsistent with 
a flat, matter-dominated universe (for which W = and 
H to = 2/3) while it is satisfied by all the ill = 2 models 
in Tab. 1. 

Supernovae la data, taken at face value, point towards 
a more negative equation of state than those listed in 
Tab. 1, typically in the W ^ —0.6 range [2]. However, 
the analyses have been done assuming two fluids with 
different equations of state, i.e. matter (W = 0) and 
'quintessence' (W — W x ), whereas in the present case 
the two fluids scale with the same, negative, equation 
of state. Since quintessence begins to dominate the en- 
ergy density quite recently, the negative equation of state 
takes over later than in our model, where it has been neg- 
ative since much previous epochs. As a consequence, the 
supernovae bounds on the equation of state of the uni- 
verse should be somehow relaxed in our model compared 
to quintessence. 

Structure formation in a generic model of the type of 
eq. (21) was studied in [13] (where also the baryonic com- 
ponent was considered). There, it was shown that the 
non-zero - and constant- SIdm allows the growth of per- 
turbations even if the expansion is accelerated. This be- 
havior is to be contrasted with the usual quintessence 
or cosmological constant case, where the perturbations 
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freeze out soon after the takeover because in that case 
SJm drops quickly to zero. 

The model we have presented in this paper is just the 
simplest version of a large family. In order to discuss 
its phenomenological implications thoroughly, one should 
specify it in detail. Different dependences of the dark 
matter mass on the small radius could also be envisaged, 
possibly leading to different predictions for the equation 
of state. Also, the couplings of the KK sector to the SM 
are crucial. If the dark matter candidate is not a SM 
singlet, then its mass variation would induce variations 
in the gauge coupling constants. Imposing the existing 
strong constraints on the latter would imply that the 4>2 
field was frozen up to a not too distant past and that the 
attractor became effective quite recently. 

It is, anyway, not trivial at all that the simplest choices, 
that is, tree-level dimensional reduction and the assump- 
tion of a O(10~ 4 /rf) Casimir energy, lead to a model 
which solves the cosmic coincidence problem and is -if 
not fully compatible with- at least very close to the obser- 
vations. A non-trivial point to consider is the persistence 
of the behavior of the tree-level solution once radiative 
corrections are included. The issue can be treated con- 
sistently only once a detailed model is given, however we 
refer the reader to the discussion in ref. [7] , where super- 
symmetry in the bulk is argued to cancel all dangerous 
corrections larger than 0(l/(rf) 2 M p ) to the radion mass. 

A model-independent prediction of the framework out- 
lined in this paper is the presence of extra dimensions in 
the 10 — 100 ^m range. Whereas for the solution of the 
hierarchy problem discussed in [5] this value was allowed, 
the link with the cosmological expansion discussed in this 
paper makes it a true prediction. The expected value for 
ri may be quite close to the present bound of 200 ^m, ob- 
tained from measurements of the Newton's law at small 
distances [14]. However, one cautionary remark is due on 
this point. The predicted value is determined by global 
observables, i.e. Ho and o'o/ao- Since the radius and 
DM are strongly coupled in this model, it is plausible 
that the local value of n, inside the galactic halo, would 
be somehow different. It is difficult to estimate the size 
of this effects, which would require a dedicated numerical 
study of the halo formation in this model. However, the 
sign would plausibly go in the direction of shrinking n 
inside the halo. This is because the system would use the 
additional degree of freedom represented by the varying 
mass, to decrease the gravitational potential in overdense 
regions, thus expanding r s -and then shrinking r; - with 
respect to the average universe. 

The effect of varying mass dark particles on cosmic 
perturbations and on the structure of halos will be the 
subject of a forthcoming publication [15]. 
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